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It is assumed that U atoms in U Ge2 have a number of / electrons appropriate to give them each 
a spin s = 1 as well as one extra itinerant electron which may equally well be on one or other U 
atom. The dynamical degrees of freedom are spin-s operators of localized spins and spin- 1/2 fermi 
operators of itinerant electrons. Applying hydrostatic pressure changes the bandwidths of spin-up 
and spin-down itinerant electrons in different way, which leads to decreasing of the contribution 
of the fermions to the magnetization keeping the spin-fermion interaction unchanged. In turn the 
local spin-fermion interaction leads to ferromagnetic superconductivity. The model accounts, in 
a quantitative and natural way, for the characteristics of the coexistence of superconductivity and 
ferromagnetism in UGe2, including many of the key experimental results: metamagnetic transitions, 
quantum transition from ferromagnetism to ferromagnetic superconductivity, the position of the 
highest superconducting critical temperature etc. 



PACS numbers: 74.20.Mn, 75.50.Cc, 74.20.Rp 



U Ge2 is the first example where ferromagnetism and 
superconductivity coexist 0, 0- The superconductiv- 
ity is found experimentally only in ferromagnetic phase, 
and only in a limited pressure range (p^,, p c ). There are 
two successive quantum phase transition, from ferromag- 
netism to ferromagnetic superconductivity at p' c , and at 
higher pressure p c to paramagnetism (figla). 

As the pressure is increased there is an abrupt decrease 
of the ordered moment at p x {p' c < p x < p c ) and another 
at p c (figlb) . The ferromagnetic state below p x is referred 
to as FM2 and the high pressure ferromagnetic state as 
FM10- It has been suggested that a spin and charge 
density wave might be formed in the FM2 state, due to 
the nesting of the Fermi surface, and they are responsible 
for the transition at However, neutron diffraction 

studies have not detected any static order due to a spin 
and charge density wave. Another possibility is that the 
transition at p x is a result of a novel tuning of the Fermi 
surface topology by the magnetization^. 

The temperature dependence of the magnetization in 
UGe2 is quite different from that found in weak itiner- 
ant ferromagnets. At zero pressure, above and well away 
from p x the low temperature dependence of the magne- 
tization has the form M(T)/M(0) - [1 - (T/T c ) 3 ] 1/2 @. 
Strictly speaking, \JGei is not a weak itinerant ferro- 
magnet, and the point where ferromagnetism and super- 
conductivity disappear simultaneously is not a quantum 
critical point at all. The Curie temperature T c decreases 
while the magnetization remains unchanged. For conven- 
tional weak ferromagnets the Curie temperature scales 
with magnetization. XJGe-i differs mainly in having a 
stronger spin orbit interaction that leads to an unusually 
large magneto-crystalline anisotropy with easy magneti- 
zation axis along shortest crystallographic axis. The dif- 
ferential susceptibility has been measured, since it gives 
a measure of the spectrum of the magnetic excitations. 
The main conclusion is that the differential susceptibil- 
ity is strongly anisotropic in the high pressure FM1 and 
paramagnetic phases but weakly anisotropic in the low 
pressure FM2 phaseQ. It is plausible that increasing the 



pressure, one changes the anisotropy, which in turn shifts 
the system from itinerant behaviour to a higher pressure 
phase which is dominated by localized spins. 
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FIG. 1: (a)The p-T phase diagram of UGe 2 . T c is Curie 
temperature, T s is the superconducting transition tempera- 
ture. (b) Pressure dependence of the dimensionless magneti- 
zation per lattice site. 

So far, there are no theoretical considerations of these 
complicated phenomena. The calculations have consid- 
ered the superconductivity to appear from completely 
itinerant ferromagnetic stateBB, or have been based on 
the physics of local moments |9l Il0|. 

Motivated by the experimental findings, one assumes 
that U atoms in U Gei have a number of / electrons ap- 
propriate to give them each a spin s = 1 as well as one 
extra itinerant electron which may equally well be on one 
or other U atom. The dynamical degrees of freedom are 
spin-s spin operators Sj of localized spins and spin-1/2 
fermi operators Ci a of itinerant electrons, where i denotes 
the sites of a three dimensional lattice. The dimension- 
less magnetization M = (J-ffJ-B of the system per lattice 
site at zero temperature is M = s + m where m is the 
contribution of mobile electrons. The parameter m de- 
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pends on the microscopic parameters of the theory and 
characterizes the vacuum. If, in the vacuum state, every 
lattice site is occupied by one electron with spin up, then 
rn = 1/2, and the electrons are highly localized as in the 
uranium compounds known as " heavy-fermion systems" . 
When, in the vacuum state, some of the sites are doubly 
occupied or empty, then m < 1/2 and the electrons are 
itinerant. The system approaches the internal point (IP) 
when m — » (M = s). It corresponds to the point p x 
of the phase diagram of UGe^ (figl). 

The local spin-fermion interaction leads to an effective 
four-fermion interaction which in turn leads to p-type 
magnon-induced ferromagnetic superconductivity (FM- 
superconductivity) ■ The order parameter is a spin 
anti-parallel component of a spin-1 triplet with zero spin 
projection (tJ, + \\), The transverse spin fluctuations 
are pair forming and the longitudinal ones are pair break- 
ing. The effective potential is attractive within an inter- 
val (pf — A,pf + A), around the fermi surface pf, where 
A depends on the parameters of transverse and longitu- 
dinal spin fluctuations. When the fermions contribute to 
the magnetization of the system (m ^ 0) spin-up and 
spin-down electrons have different (majority and minor- 
ity) Fermi surfaces. If the Fermi momenta pj and p^ 
lie within the interval (pf — A,pf + A) the interaction 
between spin-up electrons, which contribute to the ma- 
jority Fermi surface, and spin-down electrons, which con- 
tribute to the minority Fermi surface, is attractive. As a 
result, spin-up electrons from the majority Fermi surface 
transfer to the minority Fermi surface and form spin anti- 
parallel Cooper pairs, while spin-down electrons from the 
minority Fermi surface transfer to the majority one and 
form spin anti-parallel Cooper pairs too. The domain be- 
tween the Fermi surfaces determines the fermions' contri- 
bution to the magnetization m, but it is cut out from the 
domain of integration in the gap equation. When the 
electrons' contribution to the magnetization increases, 
the domain of integration in the gap equation decreases, 
and for some value of m = mf, respectively at magneti- 
zation M ' = s + m! the system undergoes a transition 
from FM-superconductivity to ferromagnetism (p' c point 
on figla). At IP (m — 0) the domain of integration in the 
gap equation is largest. As a result the superconducting 
critical temperature is highest when the system is at IP. 

The anisotropy modifies the spin-wave excitations 
adding a gap in the magnon spectrum. Increasing the 
gap, the pair formation as a result of magnons' exchange 
is suppressed, which in turn leads to decreasing of the su- 
perconducting critical temperature. Hence, the most ap- 
propriate assumption, which closely matches the experi- 
mental result, is that increasing the hydrostatic pressure 
one increases the magneto-crystalline anisotropy. For 
pressures above p x the contribution of the itinerant elec- 
tron to the magnetization is zero (m = 0), and the mag- 
netization is due to magnetization of the localized spins 
M = s. Hence, the transition to the paramagnetism 
M = undergoes with jump only. 

The important point in the spin-fermion theory is the 



mechanism of driving the system from a state with mag- 
netization M = s + m>sto the internal point (IP) 
(M = s). The subtle point is the spin-fermion interac- 
tion. It splits the spin-up and spin-down Fermi surfaces 
and leads to a nonzero contribution of itinerant electrons 
to the magnetization. Driving the system to the internal 
point (IP) one has to compensate this overall shift in the 
relative position of the energy bands keeping the spin- 
fermion interaction unchanged. In this paper a mecha- 
nism of compensation by means of different changes in 
bandwidths of spin-up and spin-down electrons is con- 
sidered. The Hamiltonian of the spin-fermion model is 

H = - J • Sj - J' Yl $** &i -JiJ2$i- Si - 

<i,j> <i>j> i 

* E ( d ^ d i° + h - c + U E "^"^ ~^E^ + 

<i,j>,<j i i 

F E (cppiicn+h.c.) (1) 

<i,j> 

Here cf a and Ci a are creation and annihilation operators 
for itinerant electrons, fii a = c^.c icr are density opera- 
tors, Si = 1/2^2 c~ljT aa / c ia , , where f denotes the vec- 

aa' 

tor of Pauli matrices, are the spin operators of itinerant 
electrons, and Si are spin-s operators of localized spins. 
The sums are over all sites of a three-dimensional lattice, 
< i,j > denotes the sum over the nearest neighbors, and 
p is the chemical potential. In the J-term corre- 
sponds to a direct Heisenberg exchange of localized spins 
which is ferromagnetic (J > 0). The magnitude of the 
magnetocrystalline anisotropy is given by J'. Here I fo- 
cus on uniaxial anisotropy, J' > 0, with the easy axis of 
magnetization along the z axis. The local spin-fermion 
interaction is ferromagnetic, too (J/ > 0), and the F- 
term describes the hopping of local pairs consisting of 
spin- up and spin-down electrons (F > 01 [Tl| . 

Still, the question remains whether the off-diagonal 
hopping parameters in the Hamiltonian, which involve 
orbital overlaps between neighbouring sites, would be suf- 
ficiently large in view of the fact that U — f orbitals are 
well localized. One expects that hybridization between 
Ge electrons and U — f electrons which gives an itinerant 
character to the f-electrons, leads to larger overlaps than 
pure / orbitals. 

I introduce Schwinger representation for the local- 
ized spin operators S, = 1/2^ lta T oa' h.a 1 ■, where the 

aa' ' 

bose operators satisfy the condition j ia ji. (J — 2s. The 
partition function can be written as a path integral 
over the complex functions of the Matsubara time t, 
fia i T ) > fia ( T ) an d Grassmann functions cf a (r) and Ci a (r) 
replacing the operators in the Hamiltonian Eqs.Q with 
the functions [12]]. In terms of Schwinger bosons the the- 
ory is U(l) gauge invariant, where the bose fields have 
a charge 1, with respect to gauge transformations, while 
the fermi fields are gauge invariants. 

It is convenient to introduce two spin-singlet fermi 
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fields 



*f(r) 



/2s 



/2s 



[/a(r)c i2 (T) - /i 2 (r)cii(T)] 

/i(T)Ci«r(r) 



(2) 



which are gauge variant with charge 1 and -1 with respect 
to gauge transformations. Equations J2J can be regarded 
as a SU(2) transformation[13j and the Fermi measure is 
invariant under the change of variables. An important 
advantage is the fact that in terms of the spin-singlet 
Fermi fields the spin-fermion interaction is diagonalized 
J2 Si • s, = s/2£[#f + *f - and one accounts 



for it exactly. The total spin of the system S* ot = Sj + Sj 



can be rewritten in the form 



<s* — _ 



2 

\b A +\li B T 



s, 



1 



(3) 



where Si is the spin vector of localized spins (S 2 = s 2 ), 
and Tj and are complex vectors which depend on 
Schwinger's bosons. They are orthogonal to the spin vec- 
tor Si • Ti = Si • T+ = and satisfy T 2 = T+ 2 = 0, T { ■ 
Tf = 2. The gauge invariance imposes the conditions 
< fyf+tyf >=< igf+tyf >= 0. As a result, the dimen- 
sionless magnetization per lattice site M =< (Sl ot ) z > 
reads 



M = 



< 



> 



<S?> (4) 



At zero temperature < Sf >= s and M = s 



where 



m = 1/2 < (*f+*f - ) > is the contribution 

of the itinerant electrons. 

Rewriting the Hamiltonian in terms of A and B fields, 
one obtains the following representations for Hubbard 
and pair-hopping terms 

E^^ = -^E(*f +$ f-^) 2 . 

i i 

E = £ Vf+VfVf+Vf. (5) 



<«,j> 



<«,j> 



One can decouple these terms by means of the Hubbard- 
Stratanovich transformation, introducing a real field 
TOi(r) associated with the composite field l/2(^ i — 
and complex fields u^(r) associated with 
^fl+^i^ ^ere i? stands for A or B. Then, the ac- 
tion is quadratic with respect to the fermions and one 
can integrate them out. The obtained free energy is a 
function of the composite fields and the integral over 
them can be performed approximately by means of the 
steepest descend method. To this end one sets the first 
derivatives of the free energy with respect to composite 
fields equal to zero. These are the mean-field equations. 



The solutions of the mean-field equations are assumed to 
be constants independent of the lattice sites and bonds 
m-(r) = 771, Ujj (t ) = u R , where m is the itinerant elec- 
tron contribution to the magnetization(see Eq.Q).The 
equations for m, u and the number of itinerant electrons 
n are 



„/ = - I deN(e)(f[e B (e)]-f[e A (e)}) 



2 

J deN(e)ef[e R (e)], R — A,B (6) 



D 



deN(e) (/ [e B (e)]+f[e A (e)}) 



where / [e fi (e)] is the Fermi function, N(e) is the density 

of state for band energy = — t ^ elkS an d bandwidth 

s 

D = 2zt with S a vector connecting a site to its nearest 
neighbors and z the number of nearest neighbors. In 
equations (|5J the fermion dispersions are 

e B {t k ) = (l-ju A y k -2mU- S -^-». (7) 

I assume for simplicity a flat density of states: N(e) = 
l/D, -D/2 < e < D/2. Unlike in the Stoner model, the 
model with pair-hopping term does not de pen d strongly 
on energy variation of the density of states [Ty. Now the 
system can be analytically solved at zero temperature. 

A solution with m = exists if u A and u B are nonzero 
and have opposite signs, which in turn requires g = F/t > 
4. Then the equation for the contribution of itinerant 
electrons to the magnetization m is 



2U 
Dg~ 



{n-lf 



- I 1 "' 



In -1[ 
2.9 



sJi 
2gD 



= 0. 
(8) 

The equation (JSJ has a solution m = if D = D x , 
where D x = sJi/\n— 1|. The Coulomb parameter U is 
large parameter in the theory, so one can choose it to 
satisfy 2U/Dg > (n - l) 2 /4 + 1/4. Then, the equation 
(JSJl has only one real solution. 

To match the experimental results it is most ade- 
quate to keep the parameters of the local interactions 
U and Ji, and the number of the itinerant electrons n 
fixed. I assume that hydrostatic pressure increases the 
pair-hopping at the expense of the single-electron hop- 
ping. This means, that the pair-hopping parameter F 
increases, while the hopping parameter t decreases when 
the pressure increases. At pressure p — p x D = D x 
and F = F x , where F x = AU/z({n - l) 2 + 1). The first 
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condition is necessary to have a zero to solution, the sec- 
ond one ensures an abrupt decrease of magnetization at 
p x . The last assumption is that when the pressure in- 
creases the parameter F scales like 1/D, more exactly 
F/F x — D x /D. Above ^ the parameters of the itinerant 
electrons remain unchanged. One can find justification 
of this assumption in the experimental fact that above 
p x the physics of the system is dominated by the local- 
ized spins. It is important to stress that the transition, 
to the paramagnetism M = 0, undergoes with jump, be- 
cause above p x the magnetization results from the local- 
ized spins. The contribution of the itinerant electrons to 
the magnetization m as a function of the pair-hopping 
parameter F/F x is depicted in fig. 2 for (n — l) 2 = 0.2, 
and 1.6U — 3sJ;. 
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FIG. 2: The contribution of the itinerant electrons to the 
magnetization m as a function of the pair-hopping parameter 
F/F x . 

The graph (fig. 2) as well as the expressions for F x and 
D x are an artifact of the approximate treatment of the 
density of states iV(e). More accurate account for the en- 
ergy dependence of N (e) will give us different conditions 
for the parameters and more realistic dependence of the 
magnetization on the parameters. 

The proposed model of U Ge2 differs from the models 



discussed in 

011 IS 

IToj in many aspects. First, degrees 
of freedom associated with localized spins and itinerant 
electrons are introduced, which enables one to describe 
two different ferromagnetic phases FM1 and FM2 fig. 1. 
The resistivity measurements reveal the presence of an 
additional phase line that lies entirely within the ferro- 
magnetic phase. It is suggested by a strong anomaly 
seen in the resistivity pi Il4| . The characteristic temper- 
ature of this transition, T x (p), decreases with pressure 
and disappears at a pressure p x (IP) at which the su- 
perconductivity is strongest. For pressures below p x the 
itinerant electrons contribute to the magnetization, while 
for pressures above p x the ferromagnetism is dominated 
by localized spins. This suggests to define T x by the 
equation m(T x ) = 0. Above T x the itinerant electrons 
do not contribute to the magnetization, and the ferro- 
magnetism is entirely dominated by the spin fluctuations 
of the localized spins, while below T x the itinerant elec- 
trons take part in the formation of the spin fluctuations. 
In particular, the itinerant electron mass renormalization 
is different below T x and above this temperature. As a re- 
sult, the slope in the dp/dT versus T diagram is different 
above T x and well below T x . Increasing the tempera- 
ture from below T x the slope changes smoothly from its 
value well below T x to its value above T x , This means 
a non-Fermi-liquid temperature dependence of the resis- 
tivity within a temperature interval around T x 14] . The 
present description of the ferromagnetism above and be- 
low T x is in very good agreement with the experimental 
finding that the high pressure ferromagnetic phase might 
have the more localized character 0. 

Second, the model explains in a unified way the su- 
perconductivity and T x transition near p x point. At p x 
the contribution of the itinerant electron to the magne- 
tization to becomes equal to zero and hence it is the end 
of the T x line, as follows from the definition above. On 
the other side, it was explained that when m = the 
superconducting critical temperature is highest. 
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